Elliptic flow of thermal dileptons as a probe of QCD matter 
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We study the variation of eUiptic flow of thermal dileptons with transverse momentum (pt) and 
invariant mass (M) for Pb+Pb collisions at ^/Snn = 2.76 TeV for 30-40% centrality. The dilepton 
productions from quark gluon plasma (QGP) and hot hadrons have been considered including the 
spectral change of light vector mesons in the thermal bath. The space time evolution has been 
carried out within the framework of 2+1 dimensional ideal hydrodynamics with lattice+hadron 
resonance gas equation of state. We find that a judicious selection of AI window can be used to 
extract the collective properties of quark matter, hadronic matter and also get a distinct signature of 
medium effects on vector mesons. Results within the ambit of the present work for nuclear collisions 
at Large Hadron Collider (LHC) indicate a reduction of elliptic flow («2) for M beyond (f> mass, 
which if observed experimentally would imitate the small momentum anisotropy in the early stage 
of the collective motion of the partons. We also observe that the magnitude of the elliptic flow is 
significantly larger at LHC than at Relativistic Ifeavy Ion Collider (RHIC) collision conditions. 
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Collision between nuclei at relativistic energies 
provides an opportunity to study Quantum Chromo- 
dynamics (QCD) at non-zero temperatures and den- 
sities. Calculations based on lattice QCD (LQCD) 
predict that at high temperatures and/or densities 
hadronic matter melts down to a state of decon- 
fined thermal quarks and gluons - a phase of mat- 
ter called quark gluon plasma (QGP). The weakly 
interacting picture of the QGP stems from the per- 
ception of asymptotic freedom of QCD at high tem- 
peratures and densities. However, the LQCD cal- 
culations indicate that the deconfined system does 
not reach the Stefan-Boltzmann limit corresponding 
to the non-interacting massless gas of quarks, anti- 
quarks and gluons even at a temperature as high as 
1000 MeV [H. The RHIC and LHC are the two ex- 
perimental facilities where the QGP can be created 
by colliding nuclei at high energies. Several probes 
to study the properties of QGP have been proposed. 
Among those the azimuthal anisotropy or the elliptic 
flow {V2) of the produced particles have been shown 
to be sensitive to the initial condition and the equa- 
tion of state (EoS) of the evolving matter formed in 
heavy ion collision (HIC)(see for review). 

Non-central heavy-ion collisions provide an 
anisotropic spatial configuration which together 
with the interactions among the constituents develop 
pressure gradients of different magnitude along dif- 
ferent spatial directions. With the expansion the 
spatial anisotropy reduces and the momentum space 
anisotropy builds up rapidly. The V2 is a measure of 
this momentum space anisotropy which is defined 
as: V2 = (cos(2(0 - *))) = {pI - P^) / {pI + pD , 
where Px and Py are the x and y component of the 
particle momenta, (j) is the azimuthal angle of the 



produced particles and is the angle subtended by 
the plane containing the beam axis and impact pa- 
rameter with x-direction. Comparison of measured 
V2 with those calculated using relativistic hydrody- 
namic and transport approaches have lead to several 
important results. The most important of these is 
the small shear viscosity to entropy ratio of the QGP 
compared to other known fluids [5|] . The mass order- 
ing of V2 of identified hadrons, clustering of V2 sepa- 
rately for baryons and mesons at intermediate pt are 
considered as signatures of partonic coalescence as a 
mechanism of hadron production @, 01 ■ In contrast 
to hadrons, which are predominantly emitted from 
the freeze-out surface of fireball, the electromagnet- 
ically interacting particles (real photons and lepton 
pairs) are considered as penetrating probes fsl (see 
Refs. (9I4T1I for review) which can carry informa- 
tion from the hot interior of the system. Therefore, 
the analysis of V2 of lepton pairs and photons can 
provide information of the pristine stage of the mat- 
ter produced in HIC. The V2 of dileptons [12] and 
photons [m have been evaluated for RHIC energies 
and shown that it can be used as effective probes to 
extract the properties of the partonic plasma. The 
sensitivity of V2 of lepton pairs on EoS has been elab- 
orated in [1^ for RHIC collision conditions. The lep- 
ton pairs are produced from each space time point 
of the system and hence the study of V2 of lepton 
pairs will shed light on the time evolution of col- 
lectivity in the system. Further, in addition to px, 
dileptons have an additional kinematic variable, M 
unlike photon. The evolution of radial fiow can be 
estimated by studying the pt spectra of lepton pairs 



for different M windows 15, [Tg. The radial flow al- 
ters the shape of the px spectra of dileptons - it kicks 
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the low pt pairs to the higher px domain, making 
the spectra flatter. Therefore, the presence of large 
radial flow may diminish the magnitude of V2 at low 
Pt [i3 and this effect will be larger when the ra- 
dial flow is large i.e. in the hadronic phase which 
corresponds to lepton pairs with M ~ nip. 

It has been argued that the anisotropic momen- 
tum distribution of the hadrons can bring the in- 
formation on the interaction of the dense phase of 
the system [2j despite the fact that the hadrons are 
emitted from the freeze-out surfaces when the sys- 
tem is too dilute to support collectivity. Therefore, a 
suitable dynamical model is required to extrapolate 
the final hadronic spectra backward in time to get 
the information about the early dense phase. Such 
an extrapolation is not required for lepton pairs be- 
cause they are emitted from the entire space-time 
volume of the system. Therefore, the V2 of lepton 
pairs provide information of the hot and dense phase 
directly. The V2 of dileptons can also be used to 
to test the validity and efficiency of the extrapola- 
tion required for hadronic V2 ■ We will also see below 
that the pT integrated M distribution of lepton pairs 
with M (> m^) originate from the early time, pro- 
viding information of the partonic phase and pairs 
with M < nip are chiefly produced later from the 
hadronic phase. Therefore, the pr integrated M dis- 
tribution of lepton pairs may be used as a chronome- 
ter of the heavy ion collisions. On the other hand, 
the variation of V2 with px for different M windows 
may be used as a flow-meter. 

In this work we study the M and pr dependence 
of V2 for dileptons using relativistic ideal hydrody- 
namics [Tsj assuming boost invariance along the lon- 
gitudinal direction for LHC collision conditions. 
We argue that the V2 of lepton pairs from QGP and 
hadronic phases can be estimated with an appropri- 
ate choice of M. In addition to this we also show 
that the thermal effects on the p spectral function 
is visible through V2 in the low mass (below p peak) 
dileptons. 

The elliptic flow of dilepton, V2{pt, M) can be de- 
fined as: 

V2 = J ^ T (1) 

Z)/ \d^pTdM2dy\y='^) '^^ 

where the ^ stands for summation over Quark Mat- 
ter(QM) and Hadronic Matter(HM) phases. The 
quantity dN j (PprdM'^ dy\y=o appearing in Eq.[T]can 
be obtained from the dilepton production per unit 
four volume, dN / d'^pd'^x in a thermalized medium 
by integrating over the space-time evolution of the 
system. The quantity, dN / d^pd^x is given by 

(2) 



where W^^{po,p) is the electromagnetic (EM) 
current correlator, g^'^ is the metric tensor, a 
is the electromagnetic coupling, p = ipo^P) 
is the four-momentum of the pair, T is the 
temperature of the thermal bath, L{M'^) = 

arises from the Dirac 
spinors (lepton pair) in the final state ,M is the in- 
variant mass of the lepton pair and mi is the lepton 
mass. 

For QGP Eq. [2] leads to the standard rate of lep- 
ton pair production (in the leading order) from an- 
nihilation of qq pairs For the dilepton produc- 
tion from hot hadrons we briefly outline the dilepton 
production processes here and refer to f2T'| for de- 
tails. For the low mass dilepton production from HM 
the decays of thermal light vector mesons namely 
p, uj and (j) have been considered. The change of 
spectral function o f p due to its interaction with 
TT, w, ai , /ii (see [2l|, [2§| for details) and baryons [ll] 
have been included in evaluating the production of 
lepton pairs from HM. For the w spectral function 
the width at non-zero temperature is taken from 
Ref [l^l and medium effects on cj) is ignored here. 
The continuum part of the spectral function of p 
and UJ have also been included in dilepton produc- 
tion rate [l^, H^- In the present work dileptons 
from non-thermal sources e.g. from the Drell-Yan 
process and decays of heavy flavours |2^ have been 
ignored in evaluating the elliptic flow of lepton pairs 
from QM and HM. If the charm and bottom quarks 
do not thermalize then they are not part of the flow- 
ing QGP and hence do not contribute to the elliptic 
flow. The model employed in the present work leads 
to a good agreement with NA60 dilepton data [13] 
for SPS collision conditions [28| . 

To evaluate V2 from Eq. [T] one needs to integrate 
the production rate given by Eq. [2] over the space 
time evolution of the system - from the initial QGP 
phase to the final hadronic freeze-out state through 
a phase transition in the intermediate stage. We as- 
sume that the matter is formed in QGP phase with 
negligible net baryon density. The initial condition 
required to solve the hydrodynamic equations for 
the description of the matter produced in Pb+Pb 
colHsion at y^s^ = 2.76 TeV for 30-40% central- 
ity are taken as follows: Ti ~ 456 MeV is the 
value of the initial temperature corresponding to 
the maximum of the initial energy density profile 
at the thermalization time Ti = 0.6 fm/c. The value 
of the transition temperature, Tc for quark hadron 
conversion is taken as 175 MeV. The EoS required 
to close the hydrodynamic equations is constructed 
by complementing Wuppertal-Budapest lattice sim- 
ulation [l[ with a hadron resonance gas compris- 
ing of all the hadronic resonances up to mass of 
2.5 GeV [H US]. The energy of the lepton pair 
(po) originating from a hydrodynamically expand- 
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ing system should be replaced by its value {p ■ u) 
in the co-moving frame which is given by: p ■ u = 
7t(Mt cosh(?/ — 77) — VxPt cos (j) — VyPr sin (/>), where 
u = (7, 'fv), is the fluid four- velocity, y is the rapidity 



defined as: 



The variations of 



-1/2 



and Vy are the x and y component 



and rj is the space-time rapidity, 7t = (1 — 



of the velocity. The system is assumed to get out of 
chemical equilibrium at T = T^h = 170 MeV [U. 
The kinetic freeze-out temperature Tp = 130 MeV 
is fixed from the pr spectra of the produced hadrons 
at the same collision energy of Pb-|-Pb system. The 
EoS and the values of the parameters mentioned 
above are constrained by the pt spectra (for — 5% 
centrality) and elliptic flow (for 10 — 50% centrality) 
of charged hadrons [2^ measured by ALICE collab- 
oration [s^l ■ 
% 




FIG. 1: Left panel: Variations of spatial and momen- 
tum space anisotropy with proper time. Right panel: 
Constant temperature contours denoting the space time 
boundaries of the QGP and hadronic phases. 
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FIG. 2: Left panel: Invariant mass distribution of lep- 
ton pairs from quark matter and hadronic matter (with 
and without medium effects). Right panel: Fractional 
contribution of lepton pairs for various invariant mass 
windows as a function of average proper time (see text 
for details). 



The spatial anisotropy of the matter produced in 
non-central HIC is defined as: ex = ^pr^^j^, where 
the average is taken over the energy distribution of 
the system formed in a collision of given impact pa- 
rameter. The resulting momentum anisotropy of the 
interacting system produced in non-central HIC is 



~ J dx J dy [Tyv+T'"'] • 

(e^^) and (cp) with time have been displayed in the 
left panel of Fig. [U for 30-40% centrality. We find 
that Cx (ep) reduces (increases) with time due to the 
prevailing pressure gradients which causes the mat- 
ter to flow. Since V2 oc Cp therefore, the variations of 
f 2 and e with r are expected to be similar. In Fig. [T] 
(right) we depict the constant temperature contours 
corresponding to Tc = 175 MeV and T/ = 130 MeV 
in the t — x plane (at zero abscissa) indicating the 
boundaries for the QM and HM phases respectively. 
The life time of the QM phase ^ 6 fm/c and the 
duration of the HM is~ 6—12 fm/c. Throughout 
this work by early and late will approximately mean 
the duration of the QM and HM respectively. 

With all the ingredients mentioned above we eval- 
uate the Pt integrated M distribution of lepton pairs 
originating from QM and HM (with and without 
medium effects on the spectral functions of p and lu). 
The results are displayed in Fig. [5] (left) for the initial 
conditions and centrality mentioned above. We ob- 
serve that for M > the QM contributions dom- 
inate. For Mp < M < the HM shines brighter 
than QM. For M < Mp, the HM (sohd line) over 
shines the QM due to the enhanced contributions 
primarily from the medium induced broadening of p 
spectral function. However, the contributions from 
QM and HM become comparable in this region of 
M if the medium effects on p spectral function is ig- 
nored (dotted line). Therefore, the results depicted 
in Fig. [2] (left) indicate that a suitable choice of M 
window will enable us to unravel the contributions 
from a particular phase (QM or HM). 

To further quantify these issues we evaluate the 
following quantity: 



F = 



dN 



d^xd^PTdM-dy ) dxdydT^Tdrd^pTdAP 



f ( <iN_ 

J yd'^xd^pTdl 



dM^a 



dxdydrjTdrd^ prdM'^ 



(3) 



where the M integration in both the numerator and 
denominator are performed for selective M windows 
from Ml to A/2 with mean M defined as (M) = 
{Ml + M2)/2. The prime in J' in the numerator 
indicates that the r integration in the numerator is 
done from ri — Ti to T2 — Ti + Ar with progres- 
sive increment of At, while in the denominator the 
integration is done over the entire lifetime of the 
system. In the right panel of Fig. [2l is plotted 
against t^v{= {ti +T2)/2). The results substantiate 
the fact that pairs with high {M) ^2.5 GeV origi- 
nate from QM (rat, ^ 6 fm/c, QGP phase) and pairs 
with (M) ^ 0.77 GeV mostly emanate from the HM 
phase (rai, > 6 fm/c). The change in the proper- 
ties of p due to its interaction with thermal hadrons 
in the bath is also visible through F evaluated for 
{M) ~ 0.3 GeV with and without medium effects. 
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This clearly indicates that the (A/) distribution of 
lepton pairs can be exploited to extract collectivity 
of different phases of the evolving matter. 
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FIG. 3: [a] and [b] indicate elliptic flow of lepton pairs 
as a function of pr for various M windows, [c] displays 
the effect of the broadening of p spectral function on the 
elliptic flow for (M) = 300 MeV. [d] shows the variation 
of Rq (see text) withpT for (M) — 0.3 GeV (solid circle), 
0.77 GeV (line) and 2.5 GeV (open circle). AU the results 
displayed here are for 30-40% centrality. 



Fig- 131 ([a] and [b]) show the differential elliptic 
flow, V2{pt) of dileptons arising from various {M) 
domains. We observe that for (M) = 2.5 GeV, V2 
is small for the entire pr range because these pairs 
arise dominantly from the QM epoch (see Fig. [21 
right panel) when the flow is not developed fully. 
By the time (6-12 fm/c) when the pairs are emitted 
predominantly in the region (M) = 0.77 GeV, the 
flow is fully developed which gives rise to large V2 ■ It 
is also interesting to note that the medium induced 
enhancement of p spectral function provides a vis- 
ible modification in V2 for dileptons below p peak 
(Fig. [3] [c]). The medium- induced effects lead to an 
enhancement of V2 of lepton pairs which is culmi- 
nating from the 'extra' interaction (absent when a 
vacuum p is considered) of the p with other ther- 
mal hadrons in the bath. We note that the differ- 
ential elliptic flow, V2{pt) obtained here at LHG is 
larger than the values obtained at RHIC [H, [l3| for 
all the invariant mass windows. In Fig. [3] [d] we 
depict the variation of Rq with pT for (M) — 0.3 
GeV (line with solid circle) 0.77 GeV (solid line) 
and 2.5 GeV (line with open circle). The quantity 



Rq (Rh) is defined as, Rq 
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where v. 



QM 



and vf^^ 



are the elliptic flow of QM and HM respectively. 
The results clearly illustrate that V2 of lepton pairs 
in the large {M){= 2.5 GeV) domain (open circle 
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FIG. 4: (Color online) Variation of dilepton elliptic flow 
as function of (M) for QM, HM (with and without 
medium effects) and for the entire evolution. The sym- 
bol * indicates the value of V2 for hadrons e.g. tt, kaon, 
proton and (j). 



in Fig. [3] [d]) originate from QM for the entire pT 
range considered here. The value of Rq is large in 
this domain because of the large (negligibly small) 
contributions from QM (HM) phase. It is also clear 
that the contribution from QM phase to the elliptic 
flow for {M){~ 0.77 GeV) is very small (solid line in 
Fig.[2] [d]). The value of Rr for (M) = 0.77 GeV is 
large (not shown in the figure). 

The V2 at the HM phase (either at p or (p peak) is 
larger than its value in the QGP phase (at (M) = 2.5 
GeV, say) for the entire pr range considered here. 
Therefore, the pt integrated values of V2 should also 
retain this character at the corresponding values of 
(M) , which is clearly observed in Fig. S] which dis- 
plays the variation of V2{{M)) with (M). The V2 
(oc Ep) of QM is small because of the small pres- 
sure gradient in the QGP phase. The V2 resulting 
from hadronic phase has a peak around p pole in- 
dicating the full development of the flow in the HM 
phase. For (A/) > the V2 obtained from the 
combined phases approach the value corresponding 
to the V2 for QGP. Therefore, measurement of V2 for 
large (M) will bring information of the QGP phase 
at the earliest time of the evolution. It is important 
to note that the pt integrated V2{{M)) of lepton 
pairs with (M) ~ mTr,mK is close to the hadronic 
V2 and V2 (symbol * in Fig. S]) if the thermal effects 
on p properties are included. Exclusion of medium 
effects give lower V2 for lepton pairs compared to 
hadrons. The fact that the V2 of the (penetrating) 
lepton pairs are similar in magnitude to the V2 of 
hadrons for {{M) ^ mT^,mK, mproton etc), it ascer- 
tains that the anisotropic momentum distribution of 
hadrons carry the information of the HM phase with 
duration ^ 6—12 fni/c (left panel of Fig. [1]). We also 
observe that the variation of V2{{M)) with (A/) has 
a structure similar to dN/dM vs M. As indicated 
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by Eq.fflwe can write V2{{M)) - T,i=qm.hm ^2 x fh 
where fi is the fraction of QM or HM from various 
space-time regions. The structure of dN/dM is re- 
flected in V2{{M)) through fi. We find that the mag- 
nitude of W2((M)) at LHC is larger than its value at 
RHIC. 

In conclusion, we have evaluated the V2 of dilep- 
tons originating from the Pb-I-Pb collisions at ^ snn 
= 2.76 TeV for 30-40% centrality. Our study shows 
that V2(M) provides useful information on the col- 
lective motion of the evolving QCD matter formed in 
high energy heavy-ion collisions. If the heavy quarks 
(charm and bottom) produced in HIC do not ther- 
malize and hence not become part of the flowing 
QGP then the lepton pairs which originate from the 
decays of heavy flavours will not contribute to the V2 
of lepton pairs. In such a scenario an experimental 
observation of the reduction of V2 (M) with increas- 
ing M beyond (f) mass would reflect the presence of 



small momentum space anisotropy through modest 
collective motion in the QM phase. We observe that 
V2{{M)) of the penetrating probe (lepton pairs) for 
(M) — TOtt and rriK is similar to the hadronic and 
V2 when the medium induced change in the p spec- 
tral function is included in evaluating the dilepton 
spectra. The medium effects are large during the 
dense phase of the hadronic system, therefore, this 
validates the findings that the hadronic V2 carry the 
information of the dense part of the hadronic phase. 
Our study also establishes the fact that the invariant 
mass dependence of dilepton V2 can in principle act 
as a clock for the space time evolution of the system 
formed in HIC. 
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